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1 Introduction 

Inflation is the dominant paradigm to explain the approximate isotropy and homogeneity 
of the universe. It also gives rise to quantum perturbations which lead to the observed 
anisotropy in the Cosmic Microwave Background, and which seed the formation of large 
scale structure in the universe. During inflation, space-time is well approximated by four 
dimensional de Sitter space, which is a maximally symmetric FRW cosmology, with the 
symmetry group 50(1,4). The time evolution of the inflaton and its back-reaction on the 
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metric breaks these symmetries, but this breaking is small if the slow roll conditions are 
satisfied. 

We will refer to the symmetry of de Sitter space as the conformal symmetry 

group. It includes translations and rotations along the spatial directions, as well as a scale 
transformation, and three special conformal transformations. In this paper, we explore the 
constraints imposed by these 50(1,4) symmetries on the perturbations produced during 
inflation. More specifically, we derive Ward identities arising due to the scale and special 
conformal transformations for the correlation functions of these perturbations. Our Ward 
identities incorporate the breaking of the 50(1,4) symmetry as well, and are valid to all 
orders in the slow roll parameters. 

The analysis we carry out is based on symmetries alone, and is independent of specific 
models. As a result, the Ward identities we obtain can provide robust model independent 
checks of the central idea behind a large class of inflationary models, namely, that the 
inflationary dynamics (including the scalar sector) preserves approximate conformal invari¬ 
ance. These results should apply not only to slow roll models with different shapes of the 
inflationary potential, but also in situations where higher derivative corrections can become 
important, such as in string theory scenarios, with the Hubble scale during inflation being 
of order the string scale, which in turn is much smaller than the Planck scale. 

The 50(1,4) symmetry is also the symmetry group of a 3-dimensional Euclidean con¬ 
formal field theory, which is the motivation behind our calling it the conformal group. 
However, we should mention at the outset that we do not assume a dS/CFT type of corre¬ 
spondence in deriving our results. Rather the connection with a conformal field theory (with 
the breaking of conformal invariance also included) is only for the purpose of organizing 
our discussion of the symmetries. 

This paper is organized as follows. Section 2 contains the basic setup. The central ideas 
and key results behind the derivation of the Ward identities are then discussed in section 3. 
For our analysis, it is useful to work with the late time wave function of the universe, when 
the modes of interest have exited the horizon. Constraints imposed by symmetries on the 
coefficient functions determining the wave function are discussed in section 4. The late time 
behaviour of the modes in the canonical slow roll model of inflation is discussed in subsection 
6.1, and some aspects which arise when higher derivative corrections are incorporated are 
discussed in subsection 6.2. We end with conclusions in section 7. The three appendices 
contain important supplementary material. 

The analysis we carry out is based on the seminal works [1] and [2]. It also develops ideas 
earlier reported in [3], [4] and [5]. There are many other references also of relevance. The 
use of conformal symmetry to constrain inflationary correlation functions has also been dis¬ 
cussed in [6-19]. Approaches where the conformal symmetries are often thought of as being 
non-linearly realized include [20-39]. The idea of using time and spatial reparametrizations 
to derive Ward identities in the context of AdS was first discussed in [40] . 

Notation: Before proceeding, let us clarify the notation we will follow in this paper. 
A dot above a quantity represents a time derivative, e.g. (j) = d(p/dt. Spatial three vectors 
are written in boldface, e.g. x,k, etc. Also, etc. represent the magnitudes of the 

vectors ka,kh, whereas ki,kj, etc. represent the components of k. Unless otherwise 
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stated, the spatial indices i,j, etc. will be raised and lowered using the Kronecker delta, 
5ij. 


2 Essential Ideas 

In this section, we will outline the essential ideas behind the derivation of the Ward identi¬ 
ties. Our discussion will be general and not tied to any specific model. In sections 6.1 and 
6 . 2 , we will discuss the concrete cases of the canonical model of slow roll inflation, and the 
presence of higher derivatives, respectively. 

The dynamical degrees of freedom in the theories we consider will be the metric and a 
single scalar field We work with the ADM form of the metric, 

ds^ = -N^dt'^ + hj (dx^ + N^dt){dx^ + Wdt), (2.1) 

with N and being the lapse and shift functions respectively. We choose the gauge 

7V = l,iV* = 0. (2.2) 

This gauge is called the synchronous gauge. 

The unperturbed background FRW solution is 

ds'^ = —dt^ + a^{t) 6ij dx'^dxK (2.3) 

The Hubble parameter is given by 

H = -. (2.4) 

a 

Including the metric perturbations, denoted by 7 *^, gives 

hij = a^{t) [5ij + 7 ij]. (2.5) 

Similarly, expanding the inflaton about the background value (^(t) gives 

4> = ^{t) + 5(1). ( 2 . 6 ) 


The gauge choice, eq.(2.2), does not fix all the coordinate reparametrization invariance. 
There are two kinds of residual gauge transformations which can be carried out. These are 
spatial reparametrizations, 


X* ^ X® -1- x®(x). 

(2.7) 

under which 


hij y hij “h ^iVj ^jVi■) 

or equivalently 

Jij lij + 

( 2 . 8 ) 

(2.9) 

®The discussion can be extended to include additional scalars. However, 
predictions are not easy to make in such models. 

model independent observational 


3 




We can also perform time reparametrizations 


t —y t -\- , 

along with accompanying spatial reparametrizations of the form 

X* ^ X* + x) 


( 2 . 10 ) 


( 2 . 11 ) 


with 


under which 


/■* 1 

w\t,x) = die{x) J dt' 


dcj) ^ dcj) + (j){t) e{x), 

lij -)■ lij + 2 6ij(-] e{x) + [diWj + djWi). 


( 2 . 12 ) 

(2.13) 

(2.14) 


Using the homogeneity of the background FRW solution, we can expand the perturba¬ 
tions in a basis of modes carrying fixed comoving momenta. Let ^ be a generic perturbation. 
Then 

/ rj'ih 

P-15) 

where the comoving momentum is k. We will be interested in the behaviour of the pertur¬ 
bations at late times, when the modes of interest have left the horizon. 


k^/a^ < 


(2.16) 


Using the time reparametrization symmetry at late times, we can set 

Scf) = 0 . 


(2.17) 


In this gauge, the perturbations freeze out once they exit the horizon, i.e., they become 
time independent, since their subsequent evolution becomes dominated by a frictional term 
proportional to the Hubble parameter. The remaining gauge invariance now corresponds 
to spatial reparametrizations, eq.(2.7). The choice of gauge eq.(2.17), and the freeze out of 
modes will be discussed in greater detail for the canonical slow roll model in section 6.1, 
and in the presence of higher derivative terms in section 6.2. 

In the gauge eq.(2.I7), all the remaining perturbations arise from the metric. We can 
decompose them as 

dij + 7ij = + 7jj], (2.18) 

where 7 jj is the traceless component. ^ determines the perturbations in the trace of the 
metric. To linear order in perturbations, we see from eq.(2.18) that 

Jij = 2C6ij+%j. (2.19) 

Going beyond the linear order, we will find that the definition given in eq.(2.18) leads 
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to a simplification in our discussion of symmetries. To be more specific, it will turn out 
that the coefficient functions for the trace of the stress tensor will transform in a canonical 
way with this choice of variables. 

It will be useful to carry out our symmetry based analysis in terms of the wave function 
of the universe. This wave function is actually afunctional of the perturbations. Expanding 
at late times, when the perturbations become time independent, we get 


= exp 


i J (fx(fyC{x)({y){T{x)T{y)) 

j d?xd^yf{x)%j{y) {T{x)T^{y)) 

i J (fxd^y%j{x)%i{y) {f^^{x)f'^\y)) 

^ J d^xd^yd^zf{x)C,{y)f{z){T{x)T{y)T{z)) 
i J d^xd^yd^zC{x)C{y)%{z) {T{x)T{y)T^ [z)) 

i J d^xd^yd^zC,{x)%(y)%i[z){T{x)f^^[y)f'^\z)) 

^ J d^xd^yd^z%j{x)%i{y)jrnniz) {f’'^{x)f’^fy)f"^'^{z)) 


( 2 . 20 ) 


771 ! n\ j ^ Xi ‘ ‘ ' d X-yfi+n C(^l) ’ ’ ’ C(^TTi) Tziji ■ ■ ■ ^injn ) X 

(r(a;i) • • • (xm+i) • • • (a;^+„)) + • • • . 

The quadratic terms in f and 7 ^ correspond to a Gaussian wave function; higher order 
terms give rise to non-Gaussianity. 

Invariance with respect to the residual gauge invariance, namely with respect to the spa¬ 
tial reparametrization eq.(2.7), imposes constraints on the coefficient functions {T{x)T[y)), 
(x)T^\y)) etc, which appear in this expansion. In fact, these coefficient functions have 
been written in a suggestive manner because the constraints take the form of Ward identi¬ 
ties which are satisfied by correlation functions of the stress-energy tensor in a conformal 
field theory. This will be discussed further in section 4. Note that in eq.(2.20) we have also 
included a mixed term between C, and for generality, although such a term will vanish on 
further gauge fixing the spatial reparametrization invariance suitably, as we will see later. 
Let us also mention that as per our conventions, eq(2.19), T is related to the trace of the 
stress tensor T^- by 

T = 2Tii = 2T, ( 2 . 21 ) 

so that the coefhcient function for a general metric perturbation, 'jij, is T*-^. Also, Tij is 
the traceless part of the stress-energy tensor T^j. 

The invariance with respect to spatial reparametrizations eq.(2.7) arises as follows. 
The wave function as a functional of the late time value for a generic perturbation ^ can 
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be written as a path integral 


( 2 . 22 ) 


m= 

J initial 


initial 


where the initial conditions will be taken to be the Bunch-Davies vacuum. The action S 
has a pre-factor 1/G ~ By suitably rescaling fields in terms of the Hubble parameter 

H, we see that 



where S contains the rescaled fields which have been made dimensionless by the rescaling. 
Since no gravity waves have been detected so far, we know that ^ 



(2.24) 


Thus the path integral on the RHS of eq.(2.22) can be evaluated in the semi-classical limit, 
by solving the equations of motion subject to the boundary conditions at late and early 
times. In particular, in the gauge eq.(2.2), the N, equations must also be imposed. These 


equations give rise to the invariance of the wave function under spatial reparametrizations, 
eq.(2.7), after hxing the gauge, eq.(2.17), at late times. 

3 Ward Identities 

We are now ready to discuss the derivation of the Ward identities. We will be interested 
in the Ward identities which arise due to scale and special conformal transformations. It is 
useful to first consider the case of de Sitter space, with the background metric 


= — dt‘^ + e^^^6ijdx^dx^. 


(3.1) 


This metric is well known to have an S'0(l,4) symmetry with ten generators. Besides the 
three spatial translations, and three rotations along the spatial directions, this symmetry 
group includes scale transformations, 


X* Ax*, t ^ t —— log(A) 


(3.2) 


and three special conformal transformations. 



(3.3) 



The scale and special conformal symmetries give rise to Ward identities on the corre¬ 
lation functions of the perturbations. In de Sitter space these identities are met exactly; 


2We take Mpi = « lO^^GeV. 
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in inflationary backgrounds, there are corrections that arise due to the evolving inflaton 
which breaks these symmetries. We will derive the resulting identities for the correlation 
functions to all orders in the slow roll parameters 


H 


(3.4) 


and 


<5 = 


H 


2HH 


(3.5) 


1 (jP' 
2lP' 


(3.6) 


The identities for the scale transformations are the analogues of the Callan-Symanzik equa¬ 
tions in field theory, which incorporate the running of the coupling constants. Similarly, 
we get identities for the special conformal transformations also incorporating the evolving 
inflaton. 

Before proceeding, let us note that in the canonical slow roll model, discussed in section 
6.1, e and ei are related as 

e = € 1 . (3.7) 


But more generally, when higher derivatives are included, they will not be related in this 
way. Also, for the slow roll conditions to hold. 


ei, (5, e <C 1. 


(3.8) 


The expectation values for the perturbations are obtained from the wave function in 
the standard manner. For example, for scalar perturbations C these are given by 

(C(®i) • • • C(®ri)) = ^ JPC] P%j] C(®i) • • • C{xn), (3.9) 

where M denotes the overall normalization factor in the path integral, 

Af = j[VC][V%]\^P (3.10) 

We will be interested in calculating these expectation values at late times, when the per¬ 
turbations of interest have frozen out. 

There is one important point which we must consider before we proceed. The sum over 
all metric perturbations pij on the RHS of eq.(3.9) is ill defined because we have not yet 
fixed the spatial reparametrization invariance symmetry. The integral on the RHS would 
diverge without fixing this symmetry. A conventional choice, which we will also make, is to 
take to be transverse, 

di% = 0, (3.11) 

besides also being traceless. With this further gauge fixing, the path integral on the RHS of 
eq.(3.9) becomes finite. Note that since pij freezes out at late times, the additional gauge 
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fixing required for eq.(3.11) can be achieved by a spatial reparametrization x® —>■ x® + e^{x) 
which preserves the synchronous gauge eq.( 2 . 2 ). 

Also note that after the additional gauge fixing, eq.(3.11), the resulting perturbations 
manifestly correspond to a scalar C, with spin 0 , and a tensor perturbation 7 *^ which has 
spin 2 , with respect to the rotations along the spatial directions. 

It is worth commenting here that the underlying reason for this further gauge fixing is 
that we are working with local correlation functions in a theory of quantum gravity. These 
correlations are well defined perturbatively about the inflationary background, but only 
after gauge fixing, as discussed above. 


3.1 Ward Identities for Scale Transformations 


Under a scale transformation 


X® ^ X® -|- Ax®, A <C 1, 

(3.12) 

and 7 jj transform as 


C —>■ C T A -|- A x®9j(^ , 

(3.13) 


% 

(3.14) 


(see appendix A). Note that the transversality condition, eq.(3.11), is preserved by this 
transformation. 

We can now consider changing variables in the path integral on the RHS of eq.(3.9), 
with C, and 77 transforming as given in eq.(3.13) and eq.(3.14), respectively. The measure 
in the path integral on the RHS of eq.(3.9) is invariant under spatial reparametrizations, 
and therefore under the change in eq.(3.12). 

Naively, on the basis of what has been discussed so far, one might also conclude that 
the wave function 4' is invariant under this transformation, leading to the condition 

('5(C(a;i)) • • • C(®n)) H- + (C(®i) • • • '5(C(®n))) = 0 , (3.15) 

where from eq.(3.13), 

(5C = A + Ax®diC. (3.16) 

This is incorrect. 

Under the transformation eq.(3.12), 7 ,^ transforms homogeneously, but Q has a ho¬ 
mogeneous and an inhomogeneous term in its transformation. The relations between the 
coefficient functions we obtain, as discussed in section 4, will ensure that terms in the wave 
function which are quadratic or higher order in the perturbations cancel amongst each other 
under this transformation. However, there is one term which arises from the leading term 
that is quadratic in C, in eq.(2.20), to begin with, which needs to be handled with care and 
does not cancel, as is also discussed at the end of section 4.2. The quadratic terms in the 
wave function include 

T~exp^-i j d^xd^yC{x)C{y) {T{x)T{y))^. 


( 3 . 17 ) 


After the transformation eq.(3.13), we get a piece arising from the inhomogeneous term in 
the transformation of C, 

C —^ C “1“ “1“ ■ ■ ■ ) (3.18) 

which will now be linear in (^, 

h'I'~exp^ —A J cfix (fy ({x) {T{x)T{y))^ . (3.19) 

This term will remain uncanceled. In contrast, the homogeneous term in the transformation 

ofC, 

C C T Xx^di^ + • • • , (3.20) 

will give rise to a term which is quadratic in this will cancel against a term coming from 
the piece of T cubic in 


Before proceeding, let us note that in eq.(2.20) there is another quadratic term, 

T~exp^-y d^xd^yC{x)%j{y) {T{x)T^{y))^, (3.21) 

involving both ^ and %j, which could also have potentially contributed an additional piece. 
However, in the gauge eq.(3.11), this term in the wave function vanishes. This follows after 
noting that symmetries require the momentum space coefficient function {k\)T{k 2 )) to 
be of the form 

(r^'(fci)T(fc2)) ~ {2T:)H\k^+k2) Qhii - (3.22) 

where /3(A:i) is a dimension 3 function of ki. 


Keeping this uncanceled linear term, eq.(3.19), gives us then the correct Ward identity 

(5(C(®l)) • • • C{Xn)) + • • • + (C(®l) • • • <5(C(®n))) 

(3.23) 


= 2Ay' d^x d^y {T{x)T{y)) {C{xj_) ■ ■ ■ C{xn)C{x)). 


We will be interested in the expectation values for with non-zero momentum. Since A is 
a constant, we can drop the piece linear in A on the LHS of eq.(3.23), leading to the Ward 
identity 


E 

a=l 


X. 


d 

dXa 


(C(®i)---C(®n)) = 2 / d^xd^y{T{x)T{y)) {C{xi)---C{xn)C{x)). (3.24) 
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Expressing this in momentum space gives ^ 

^3(71 - 1) + ^ fca (C(fci) • • • ((kr^)y = 

'' a=l 


(C(k. 


n+1 




fcn + 1—>■0 

Similarly, for correlation functions of tensor perturbations ^ij we get 

f X!/ j ■ ■ ■ linjni'^n)) = 

^ a=1 “^ 

2 J d^x d^y {T{x)T{y)) {%j^{xr) ■ ■ ■yi^j^{xr^)C{x)), 

which in momentum space takes the form 

^3(71 — 1 ) + y ^ ka {ihji (^1) ■ ■ ■ linjn (kn)) = 


fc-n+l—>■0 


(3.25) 


(3.26) 


(3.27) 


Mixed identities involving both tensor and scalar perturbations can also be similarly ob¬ 
tained. These are given by 


^3(77 1) + ^ ^ feg ^ (7*1^1 (^1) ■■■ 7imim (^m)C(^m-l-l) ••• C(^r!,)) 

1 

w {liijiikl) • • • 'yimjm{km.)C{km+l) ' ' ' C{kn+l)) 


(C(fcn+l)C(-fcn+l))' 


(3.28) 


^Tl + 1—^0 


Equations (3.25) and (3.27) are examples of Maldacena consistency conditions in the 
literature [1]. These are exact to all orders in the slow roll expansion. 

The physical picture behind these relations is easy to state. The LHS of eq.(3.24) 
is the change of the n-point correlator under an overall change of scale. Exactly such a 
transformation is generated by a scalar perturbation (("(fej^q-i) in the limit of very long 
wavelength, kn+i —t 0, leading to the identity eq.(3.25). 

Comments'. The reader will note that the scale transformation eq.(3.12) is different 
from the isometry in de Sitter space, eq.(3.2). In de Sitter space, metric perturbations and 
also perturbations for test scalars freeze out at late times and become time independent. 
Thus, in effect, the scale transformation becomes eq.(3.12). In the inflationary case, once 
we choose the gauge where eq.(2.17) is met, we cannot make any time reparametrization, 

prime ' symbol on a correlation function denotes the suppression of the overall momentum conserving 
delta function. For e.g. 

(C(fel)C(fe2)) = {2^f5\k^ + fc2) (C(fcl)C(fc2)>'. 
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eq.(2.10). Thus the only symmetries available are spatial reparametrizations. 

Similarly, the special conformal transformations, which we will consider next, eq.(3.29), 
are different from the corresponding isometries in de Sitter space, eq.(3.3). However, again 
at late times, their action on time independent fields will be the same as eq.(3.29). 

It is also worth emphasizing that our derivation of the Ward identities for scale in¬ 
variance obtained here is quite general. As mentioned above, it is valid to all orders in 
the slow roll expansion, and thus should hold even when the slow roll conditions are not 
valid. The assumptions one has used are that one can go to the gauge eq.(2.17), and that 
the remaining metric perturbations in this gauge then freeze out due to the cosmological 
expansion. The residual spatial reparametrizations are then enough to give rise to the 
Ward identities above. A similar comment will also apply to the Ward identities of special 
conformal invariance we derive next. 


3.2 Ward Identities for Special Conformal Transformations 


We next turn to the special conformal transformations, 

X* ^ X* -|- a*(x), a^(x) = —2(6 • x)x^ + (3.29) 

Here there is an important extra subtlety. Consider the transformation of C, and yjj under 
eq.(3.29) (see appendix A), 


C(x) —)■ C(®) “ 2(6 • x) + a^di({x ), (3.30) 

%j{x) %(x) + a^dmjijix) + 2Mi,n{x)jjm{x) + 2Mj^{x)%m{x), (3.31) 

where Ai^-[x) is given by 

M.\j{x) = Xihj — Xjbi- (3.32) 

It is easy to see that the transformation eq.(3.31) does not preserve the transverse gauge 
condition eq.(3.11) we have chosen for 7 .^. We must therefore carry out a compensating 
coordinate transformation 


v'‘{x) = — 


X* -h u*(a;), 

^h'^^im{x) 


32 


(3.33) 


which then restores the transversality condition on 7 *^. Under this compensating transfor¬ 
mation, Q and 7 jj transform as 


C(x) ^ C(®) 


52 


4C(®) 


2b^dajjn{x) 

32 




(3.34) 
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and 


%j{x) %{x) 


66™ 


di 


( Ijmjx) \ 

V 52 J 


+ 9j 


f %mix) y 

V 52 J\ 


66 ™ 


%k{x)dj 


(1km{x^ \ 

J 


+ ljk{x)di 


(Ikmix) \ 

v^^y. 


66™(9fc7ii(®) 


f ^km{x) \ 


+ 46™7afe(a;)5a 


f %m{x) \ 

V 92 J 


i^ij + 7ij(x)) ■ 


(3.35) 


The Ward identities then arise because of the combined transformations eq.(3.30) and 
eq.(3.34) for the transformation of C, and eq.(3.31) and eq.(3.35) for the transformation of 
7 jj. Note that the compensating transformation parameter u* itself depends on 7 *^. As a 
result, the compensating transformation becomes non-linear in the perturbations. 


Once this subtlety requiring a compensating coordinate transformation is taken care of, 
the rest of the analysis follows along similar lines to that for the scale transformation case. 
The wave function T is invariant under spatial reparametrizations, and therefore under the 
combined transformations eq.(3.29) and eq.(3.33). More correctly, this is true for all terms 
in the wave function which are quadratic or higher order in the perturbations. However, 
the inhomogeneous term in eq.(3.30), —2(6 • x), gives rise to a term in the change of the 
wave function which is linear in This term does not cancel. As a result we get a Ward 
identity for scalar perturbations of the form 


4 J d^xd^y{b-x){T{x)T{y)){C{xi)---C{xr,.)C{y)) 

+ {S^C{xi) ■ ■ ■ C{Xn)) H-h {C{xi) ■ ■ ■ S'^Cix-n)) = 0, 


(3.36) 


where denotes the complete homogeneous change in ( under eq.(3.30) and eq.(3.34), 

6^C{x) = (-2(6-a;)x^ + 6^®^) dkCix) - dkCix) - ^^%rn{x) (-3 37 ^ 

In momentum space this takes the form 


92 


(6(C(fcl)) • • • C(fcn)> + • • • + (C(fcl) • • • -5(C(fcn))> = 


-2 b- 


dk. 


d \ (C(fcl)---C(fcn+l)) 


■Tl-f-l 


J (C(fcn+i)C(-fcn+i))' 


(3.38) 




where 6{C,{k)) is given by 


6(C(fc))=4'C(fc) + 6 6™fc* 

+ 2 6™ A:* 


d^k ^ 

— ^^ak-k)^uk} 

(2^ |2 7jm{k), 


(3.39) 
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and the operator is given by 


.AVh.AU( 6 .fc)('A.^ 1 + 616 .^ 


akJ 


dk dk 


d 

dk 


Similarly, for the tensor perturbations we get the Ward identity 


(^l)) ■ ■ ■ linjni^n)) + • • • + {'Jiljl (^l) ■ ■ ■ ^{linjn ) 

d \ ( 7 * 1^1 (^l) ■■■ 7injn (^ri)C(^n4-l)) 


= -2 6 - 


dkr,+lj (C(fcr^+l)C(-^r^+l))' 

where S{%j) is the complete change in yjj in momentum space, given by 




(3.40) 


(3.41) 


S{lij{k)) =Cl%j{k) + 2M'l^{k)^jrn{k) + 2 M)^{k)%^{k) 


where 


1 


+ p {kiljmik) + kjjimik)) 


+ 66 ™ 
+ 66 ™ 

- 46™ 6, 
-46™ 


/ 


d^k' %m{k') 


(27r)3 

r d^k' k', 

J ( 2 vr)‘ 


{k'i^jkik - k') + k'j%k{k - k')) 


3 Ifc — 
d?k' k 


f d-^k' 

J ( 2 vr)' 


d^k' d?k" ika-k'^) 


^ 7j2 %'(^07«m(fc - k') 

777 1 o-bi,k ) 'yhm{k fc ) 

7 a 6 (fc") 76 m(fc - fc' - k"), 


(3.42) 


(27r)^ (27r)^ \k — k' — fc"p 
d 


MUk) = b 


7 -hi 


d 


^ dk^ dk^ 


(3.43) 


Finally, we can write the general Ward identity for the variation of a mixed correlator 
involving m tensor perturbations %j{k) and (n — m) scalar perturbations C(^) 


('^( 7*171 (^ 1 )) ■ ■ ■ limjm {km) C{km+l) • • • C{kn)) + ' ' ' 

4” (7jiji(^ 1 ) ■ ■ ■ d{^i^j^{kjfk)) C{km+i) • • • C{kn)) 

4” (7njl (^ 1 ) ■ ■ ■ li-mjm {km) d{C{km+l)) • • • C{kn)) 4- • • • 
4“ (7iijl (^ 1 ) ■ ■ ■ limjm {km) C(^m+l) ' ' ' d{C{kn))) 


(b ^ ' 

(7*171(^1) ■ ■ ■ 7 imjm{km)C{km+l) • • • C{kn+l)} 


(C(fcr.+l)C(-fcn+l))' 


(3.44) 


where 6{C{k)) is given by eq.(3.39) and 5{jij{k)) is given by eq.(3.42). 

These identities are again exact, like the ones for scale transformations derived in 
section 3.1. They are valid to all orders in the slow roll expansion, and are one of the 
key results of this paper. Note that due to the non-linear nature of the transformations 
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eq.(3.34) and eq.(3.35), the resulting identities are in fact quite complicated. We will see in 
section 5, following [4], that the identity for the four point scalar perturbations in de Sitter 
space is indeed met. 

It is important to note that both connected and disconnected contributions to the 
correlation functions may be important in the Ward identities. As mentioned in eq.(2.23), 
after suitable rescaling the action has a factor of in front of it. Since this ratio is 

large, eq.(2.24), the situation in cosmology is analogous to the large N limit in AdS/CFT, 
with playing the role of . Disconnected components in the Ward identities 

can then often dominate over connected ones. More accurately, the non-linear nature 
of the transformation in the compensating spatial reparametrization means that different 
number of helds will be present in correlation functions involved in the LHS of the Ward 
identities eq.(3.38), eq.(3.41). The suppression at large N due to additional helds can be 
compensated for by including additional disconnected components. For more details and 
an explicit example see section 5. 

Let us also note that in the Ward identity for scalar perturbations eq.(3.38), for the 
cases n = 2,3, the extra terms in S({k) due to the compensating spatial reparametrization 
can be neglected to the leading order in Ff^/Mpp The extra terms are the last two terms 
on the RHS of eq.(3.39). This will become clearer from the discussion in section 5. As a 
result, for the cases n = 2,3, the Ward identity eq.(3.38) takes the form, 



{C(fc.X(M> = -2(6.A) 


(C(fcl)C(fc2)C(fc3)) 

(C(fc3)C(-fc3))' 


fca-J-O 


(3.45) 


and 



(C(fcl)C(fc2)C(fc3)) = -2(b.^) 


(C(fcl)C(fc2)C(fc3)C(fc4)) 

(C(fc4)C(-fc4))' 


fc4—^0 


(3.46) 


For the special case 6 oc fcs in eq.(3.45) and 6 oc ^4 in eq.(3.46), we obtain the conformal 
consistency relations derived in eq.(37) of [26]. However, as we will see in section 5, for 
the case n > 4 in eq.(3.38) the extra terms in SC{k) due to the compensating spatial 
reparametrization cannot be neglected. Thus, the Ward identities obtained from eq.(3.38) 
for b oc kn-j-i , n > 4, have additional terms as compared to eq.(37) of [26]. 

This concludes our summary of some of the main results of this paper. 


4 Conditions on the Coefficient Functions 

In this section, we will show how the invariance of the wave function under the scale and 
special conformal transformations, eq. (3.12) and eq.(3.29), lead to conditions on the coef¬ 
ficient functions which are analogous to Ward identities in a conformal held theory. These 
Ward identities also incorporate the breaking of conformal invariance in the inhationary 
background. 
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We start with the Ward identities of spatial reparametrizations in general. These will 
give rise to conditions analogous to Ward identities of stress-energy conservation in a field 
theory. We then obtain the identities for scale and special conformal transformations. 

In this section, it will be convenient to work with a general metric perturbation 7 *^ 
without further imposing the transversality condition eq.(3.11). Once the constraints on 
the coefficient functions have been obtained for these general perturbations, they will lead 
to constraints on expectation values which can be calculated only after further gauge fixing, 
as explained in section 3 above. 


4.1 Spatial Reparametrizations 

Under a spatial reparametrization eq.(2.7), the perturbations ( and jij transforms as 

C ^ C + ^ ^ diVj (4.1) 

and 


'Yij t 7ij 4“ 4” 9jVi ^ daVa 4“ djV + 7j;fe dit 

2 2 \ 

4- V^dkJij - - daVa%j - - daVb^ab i^ij + %) j • 


(4.2) 


See appendix A for some details of the derivation. 

Now, for the invariance of the wave function under spatial reparametrizations, the 
terms proportional to the transformation parameter Vi that get generated because of the 
transformations eq.(4.1), eq.(4.2), must cancel with one another. Consider first such terms 
which are linear in <^. These terms are produced by the first and second terms of the wave 
function, eq.(2.20). These are 


6[-- d^xd^yC{x)C{y){T{x)T{y)) 


= -g / d-^xd-^y^diVi{x)C{y) + C{x)dyiVi{y)^{T{x)T{y)), 


(4.3) 


and 




- I d^xd^yC{x)%j{y){T(x)T^{y)) 




(4.4) 


= -2 / d^xd^yC{x)dyiVj{y) {T{x)T^{y)). 


Mutual cancellation of the terms in eq.(4.3) and eq.(4.4) produces the Ward identity^ 


d,{T\x)T{y)) + - d,{T{x)T{y)) = 0. 


(4.5) 


*Note that unless otherwise stated, di stands for the derivative with respect to x*, i.e. d/dx^. 
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Using 


(4.6) 


1 


+ -SijT, 

and eq.(2.21), this can also be written as 

di{r^{x)Tiy))=0. (4.7) 

Similarly, canceling the extra terms in the wave function which are linear in yjj gives 

di{T^^{x)T'^\y)) =0. (4.8) 


Proceeding in a similar manner, and canceling terms linear in Vi and quadratic in ^ 


gives 


di{r^{x)T{y)T{z)) = i d, [h3(ai - y)] {T{x)T{z)) + ^ dj [^^(ai - z)] (T(x)T(y)} 

-^dj(T(x)T(y)T(z)}. 


(4.9) 


Eq.(4.9) can be rewritten in terms of the complete stress-energy tensor and its trace T 
as 

di{T^{x)T{y)T{z)) = ^ dj [(5^(a; - y)] {T{x)T{z)) + ^ dj [^^(ai - z)] {T{x)T{y)). (4.10) 
Similarly, canceling terms proportional to ("Tij produces the Ward identity 
d,{T^{x)f>^^{y)T{z)) = \d,[6Hx-y)]{f>^^{x)T{z)) 


+ -dj[6Hx-z)]{T{x)T>^^{y)) 
+ ldj[6Hx-y)]{f'^^{y)T{z)) 
-Sjid,[SHx-y)]{T^Hy)T{^)), 


(4.11) 


and proportional to gives 


di{T\x)f'^\y)f'^^{z)) 


= d,[5\x-z)]{f^^{x)f^\y)) 

+ ‘^-d,[5\x-z)\{f^\y)f^^{z)) 
-25,nd,[6\x - z)\{f^\y)T^{z)). 


(4.12) 


4.2 Scale Transformations 

We now turn to deriving conditions on the coefficient functions for the invariance of the 
wave function under scale transformations, eq.(3.12). The change in Q and 7 ij under these 
transformations is given by eq.(3.13) and eq.(3.14) respectively. 

The procedure we follow is the same as outlined above. Canceling the terms linear in 
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the transformation parameter A and quadratic in C gives 


I d^z{T{x)T{y)T{z)) 


d j d 




{T{x)T{y)) + Q{T{x)T{y)), 


which in momentum space takes the form 


^hm (r(fci)T(fc2)T(fc3)) = - ( 1 (nfci)T(fc2)), 


\a=l 


(4.13) 


(4.14) 


where ka = |fca|- 

In general, the ra-point correlation function of the T operators will be related under 
scaling to the (n — l)-point correlation function through the relation 


hm (r(fci)-..r(fc^)) 




(r(fci)---r(fc„_i)). 


(4.15) 


{T^{x)T>^\y)) + iS{r\x)T^\y)), (4.16) 


kh 


Requiring the cancellation of the extra quadratic terms in 7 *^ gives us the Ward identity 

which translates in momentum space to 

d 


hm {T^{ki)T^\k 2 )T{k^)) = - 

fe3->0 


a=l 


dkn 


(r^(fci)r^'(fc 2 )). 


(4.17) 


The general form of the scaling Ward identity relating the n-point correlation function 
of (n — 1) operators and one insertion of T, to the (n — l)-point correlation function of 
the T*-^ operators is 


kn^O 


'n—1 


_a=l 


d 

dka 


(f*i^i(fci) • • • r"-i^"-i(fc„_i)). 


(4.18) 


One can also write the general scaling Ward identity relating the (re + l)-point correla¬ 
tion function involving m insertions of T*-^ and (re-|-1 — m) insertions of T, with the re-point 
correlation function of m insertions of T*-^ and (re — m) of T, as 


hm (ri^Hfci)---r*"‘'"‘(fcm)T(fc^+i)---r(fc^+i)) = 

fcn + l^O 


n 


(fci) • • • r-^-(fc^)r(fc^+i) • • • r(fc„)). 
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(4.19) 












One final comment. We began this subsection by considering terms which are quadratic 
in eq.(4.13). There is also a term which is linear in both ( and the transformation 
parameter A. Since A is spatially constant, this term has support only at zero momentum 
in the wave function and we neglect it here. However, in deriving expectation values, this 
term which is uncanceled plays a crucial role, as was discussed in section 3.1 above. 


4.3 Special Conformal Transformations 


We will now derive Ward identities for the invariance of the wave function under special con¬ 
formal transformations, eq.(3.29). The change in ^ and jij under this is given by eq.(3.30) 
and eq.(3.31). Note that, as was mentioned at the beginning of this section, we are con¬ 
sidering a general graviton perturbation here and have not fixed it to be transverse; as a 
result we do not have to worry about the fact that a special conformal transformation leads 
to the gauge eq.(3.11) not being preserved. 

We start again with terms which are linear in and quadratic in Invariance of T 
then gives 


3b-(x + y) {T{x)T{y)) - ^ [a\x)d^i + a\y)dyi] {T{x)T{y)) 

= I d^z{b-z){T{x)T{y)T{z)), 


(4.20) 


which in momentum space has the form 




{T{kr)T{k^)) 



(r(fci)r(fc2)T(fc3)) 


fcs-^O 


where is the operator 


(4.21) 


C 


b 

k 




-{b-k) 


d d \ 
dk dk) 


= {b-k) 


2 ^ A 
k dk dk'^ ) 


(4.22) 


In general, the re-point correlation function of T operators will be related to the (re —l)-point 
function under a special conformal transformation as 


1 

2 


^n-l 

E 

\a=l 


c 


ka 


(r(fci)---r(fc^_i)) 


= - b- 


d 

dkr 


(r(fci)---r(fc^)) 


kn —^0 


(4.23) 


Similarly, canceling the extra terms quadratic in yjj gives 

+ Vl] {r\x)f'^\y)) = -I d^z{b- z){r^{x)f^\y)T{z)), (4.24) 
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where the action of the operator is defined by 


r^{x) = -3 (b • rr) r^{x) + - a^ix) OmT^x) 


-Ml,{x)T^Kx)-Ml^{x)r^{x\ 


(4.25) 


with Ai\j{x) as dehned in eq.(3.32). The Ward identity eq. (4.24) can be expressed in the 
momentum space as 


[vt,+v^;^{r\k^)T^\k2)) = - [h 


d 

dk^ 


{T^{kr)f^\k2)Tiks)) 


ka^O 


(4.26) 


where the momentum space operator is dehned by 


VtT^ik) = -CiT^ik) - MUk)T^^{k) - Mtj{k)T^{k), (4.27) 

with the operator as given in eq.(4.22), and A4^j{k) as given by eq.(3.43). 

Following a similar procedure as outlined above, we get the Ward identity for special 
conformal transformations relating the re-point correlation function with (re — 1) insertions 
of r*-^ and one insertion of T, to the (re — 1) point correlation function of to be 


n—l 

E 

a=l 


V 


ka 


(ri^i(fci) • • •T*-«"-i(fc„_i)) = 



(r«i(fci) • • • r"-i^-i(/e„_i)r(fc„)) 


kn —^0 


(4.28) 


In general, the Ward identity of special conformal invariance relating the re-point cor¬ 
relation function with m insertions of and (re — m) insertions of T, to the (re -|- l)-point 
correlation function with m insertions of and (re -|- 1 — rer) insertions of T is given by 


m \ ^ y n 

+ of ^ 


a=l 

= -f 6 . 


2 

— ) 

dk^^ij 


(r«i(fci) • • • r-^-(fc^)r(fc^+i) • • • r(fc^)) 


r=m+l 

(fci) • • • r-^-(fc^)r(fc^+i) • • • r(fc„+i)) 


(4.29) 


1 + 1 - 


where the operators and are given in eq.(4.27) and eq.(4.22) respectively. 

Finally, as in the case of the scale transformations, there is one remaining term linear in 
7 jj and bj, with support at zero momentum, which is uncanceled. In calculating expectation 
values, section 3, this term will vanish once we choose the transverse gauge, eq.(3.11). 


5 Explicit Checks of the Special Conformal Ward Identities 

In this section, we present a few checks of the Ward identities of special conformal trans¬ 
formations, eq.(3.38) and eq.(3.41). Our analysis above has been to the leading order in 
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H‘^ /Mpi^ and we will verify the Ward identities to this order below. 


5.1 Basic Checks 

Consider the Ward identity eq.(3.38) for the case of re = 2. We have 


(C(fcl)C(fc2)) + <! {6 5™fei I ^l(C(fci-fc)C(fc2)7.m(fc)) 


(27r)3 


a=l 

+ 2b'^k\ f ^ {%{ki - k)jjm{k)C{k2)) ] + [ki ^ k2 


= -2 b- 


(27r)3 


(5.1) 


fcs—^0 


We will first consider the identity eq.(5.1) in the de Sitter limit. We substitute 

C = - 


(5.2) 


and take the limit 0 —)■ 0. Keeping only the leading terms, we get 


{S<p{ki)6(j){k2)) = - j ^ (<5(()(A;i - fc)<5(()(A;2) 7im(fc)) 


+ ( fci -f-)- ^2 


(5.3) 


Next, we introduce suitable factors of H/Mpi. The wave function eq.(2.20) arises from 
the action eq.(6.31), which has a factor of 1/G ~ Mp^ in front of it. Thus, after suitably 
rescaling by powers of the Hubble parameter, T will go like 


T ~ exp 




d^xd^y5(f){x)5(f){y) {0{x)0{y)) 


+ ^J d^xd^y%{x)%i{y) {T^{x)f^\y)) 

+ Y\j d^xd^yd^z64i{x)54>{y)6(j){z) {0{x)0{y)0{z)) 

+ ^/ d^xd^yd^z64>{x)64>{y)%{z) {0{x)0{y)T^{z)) ^ -| 


(5.4) 


As a result, we see that the propagators {5(f) 5(f)) or {^ij^kl) behave like /Mp^ while each 
vertex, e.g., the three point vertices {0{x)0{y)0{z)) and {0{x)0{y)T^^{z)) in eq.(5.4), go 
like Mpi/H‘^. With this, one can argue that 

{5(p5(f)%j) ^ (5.5) 

From eq.(5.5), it becomes clear that to leading order in H'^/Mp^, the correlation 
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function {5(j){k\ — k)64>{k2)%mik)) in the RHS of eq.(5.3) is suppressed compared to 
{5(f){k\)5(j){k2)) in the LHS and eq.(5.3) reduces to 




{5ct>[kr)5ct>[k2)) = 0 . 


(5.6) 


This condition is the statement of conformal invariance of the two point function {5(p{ki)6(f>{k2)) 
in de Sitter space, and it is easy to verify that it is met. 


Next, let us consider the n = 2 scalar Ward identity, eq.(5.1), to the first non-trivial 
order in the slow roll approximation (but still to the leading order in /Mpi). The terms 
proportional to (CCl) ( 77 C) in eq.(5.1) scale as and can be dropped compared 

to other terms, which scale as H'^/Mp^, in the limit H <C Mpi. Eq.(5.1) then reduces to 




mi)i;(k2)) = -2[b 


1 


a=l 

In the canonical slow-roll model we have 


d 


dk^J iCiksKi-ks)) 


- (C(fcl)C(fc2)C(fc3)) 


(5.7) 


fcs-i-O 


1 


and^ 


(C(fci)C(fc2)) = {27rfdHki + ^2) ^ ^ 


1 1 


(5.8) 


(3e - 2ri) + k^k 


a=l 


a^b a>b 


(5.9) 


Working to the leading order in the slow roll parameters e, rj as well as the scalar tilt ns, 
one finds that eq.(5.7) implies 

ns = 2(r/-3e), (5.10) 

which is the correct result, [1]. 


Note that the case n = 3 for the Ward identity eq.(3.38) in the slow roll approximation 
was discussed in detail in [5]. 


Next, we turn to the graviton two-point correlator and consider the Ward identity 
eq.(3.41). Again dropping terms which are subleading in H'^/Mpp and working in the de 
Sitter limit 0 —)■ 0, we get. 


®The result eq.(5.9) is from [1], with kt = ki + k 2 + fcs. 
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E«. 

- a=l 

= -(^27Wf^(fci)(%m(fci)%z(fc2)) +27W^^(fci)(7im(fci)%i(fc2)) 

+ 2Ai^^(fc2)(%(fci)7mz(fc2)) + 2M'l^{k2){%j{kx)%m{k2)) (5-11) 

Q}fri 

+ (fcli {^jm{ki)%i{k2)) + kij {%m{ki)%i{k2))) 

QjjTn . 

+ -j^ {k2k {%j{ki)%rn{k2)) + /C2i (%'(^1 )%m(fc2))) j • 

The two-point graviton correlator is 

{%i{k)%,{-k)Y = (5.12) 

where Pijki{k) is given in eq.(5.2) of [4], An explicit calculation then shows that eq.(5.11) 
is indeed met. Note that the last two terms on the RHS of eq.(5.11) come from the 
compensating spatial reparametrization which maintains the transverse gauge for 7 *^. 


5.2 The Scalar Four Point Function 

In this subsection, we consider the Ward identity eq.(3.38) for the case n = 4. We will work 
to the leading order in H'^/Mpi, and to the leading order in ^/H, i.e. in the de Sitter limit. 
Eq.(3.38) for the case of n = 4 then gives 


a=l 


(C(«!i)c(*!2)c(fc3)c(fc4)) + q 26'“fc;y^px 


3(C(fci - k)C{k2)C{k3)C{k4)%m{k)) + {%j{ki - k)ljm{k)C,{k2)C{k3)C,{k4)) 


T (^1 ^2) T (^1 ^ ^3) T (^1 ^4) 


(5.13) 


= -2 b 


d \ 


dkf 




fes—)-0 


We next write eq.(5.13) in terms of 64> by using eq.(5.2), and take the de Sitter limit (/>—>■ 0. 
In this limit, the terms in eq.(5.13) that survive are 




{5ct>{k^)6ct>{k2)6ct>{k3)6ct>{k^)) 


T (^1 ^2) T (^1 ^3) T (^1 ^ k^j >. 


(5.14) 
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Introducing suitable factors of H/Mpi by rescaling the wave function, eq.(5.4), we see 
that for connected correlators, 

{5(t)5(t)5(t)64>) ~ ~ (5.15) 

From eq.(5.15), it seems that for the Hubble scale being much small compared to the 
Planck scale, H <C Mpi, the correlation function — k)6(j){k2)d(j){ks)5(f){k4)^im{k)) 

in the RHS of eq.(5.14) is suppressed compared to {5(p{k4)5(p{k2)64>{ks)6(f>{k4)) in the 
LHS. However, one should also consider disconnected contributions to the RHS of eq.(5.14) 
which may contribute to the same order of H/Mpi as the LHS. In particular, there is a 
disconnected contribution to the five point function {5(j)5(j)5(p64rfij) in eq.(5.14), which goes 
as 

{5(t)5(t)) {5(t)54>%j) ~ (5.16) 

and which is of the same order as {64>5(f)6(l)5(p). With these considerations, eq.(5.14) in the 
limit H <C Mpi becomes 




{6(t>{k4)6(t>{k2)64>{k^)64>{k4)) 


+ (fc2 ^ 3 ) 4“ (^2 ^ k^j I + (fci o fc2) 4“ (^1 ^ 3 ) 4“ (fci ^ ^ 4 ) 


(5.17) 


It is important to note that there are other possible disconnected contributions to the 
five point correlator {5(j){k4 — k)5(p{k2)S4>{ks)S4>{k4)jimik)), such as 


{5(p{k2)5(p{k3)){%jn{k)5(p{k-i_ -k)5(p{k4)). (5.18) 


However, this requires ^2 4- ^3 = fci 4- ^4 = 0, and will not contribute in general. 

Eq.(5.17) gives the change in the four point correlator {6(j)5(j)64>5(j)) under a special 
conformal transformation in the exact de Sitter limit. Using the relations 


and 


{5(^{k4)5(^{k2)) 


{27ry5^{k4 + k2){6^{k4)6^{k2)y 


{27rf6^{ki + k2 


1 

WiW 


(5.19) 


{lim{k)5(p{k3)6(p{k4)) = - 2 {54){k3)6(t){-k3)y{64>{k4)6(j){-k4)y{jim{k)%i{-k)y 

{ni{k)0{k3)0{k4)) 

{%rn{k)%ii-k)y{fki{k)0{k3)0ik4)) 

Mj,, {2kl){2kl) 
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we can write eq.(5.17) as 




{d(j){k^)6^{k2)6^{ks)d(P{k4)) = 


6b^Pi^ki{ki + k2) i fku k2^\ {fkl{kl + k2)0{ks)0{k4)) 

M% |fci + fc2p \ykl kf) { 2 kl){ 2 kl) 

k^ {fkijks + k4)0{ki)0{k2)) 

\kl k^J {2kl){2kl) 


+ 


k2 o k^ 


+ 


^2 ^ ^4 , 


(5.21) 


where we have used the eq.(5.12). 

Eq.(5.21) gives us the change in the scalar four point function under a special conformal 
transformation. We can verify this by performing an explicit check. The four point function 
{6(j)6(j)5(p64>) was calculated in [4]. It is given by 


{dc/){k^)ScP{k2)S(P{ks)6(^{k4)) = j W]<5</.(fci),5(/.(fc2)50(fc3)50(fc4)P[50], 


(5.22) 


where T’[(5</>] is the probability distribution function 

P[5(j)] = J[V%] . 


(5.23) 


An explicit expression for P[6(p] was obtained starting from the wave function, eq.(5.4), in 
eq.(5.3) of [4]. It is given by® 


P[d(p] = exp 


^'kpi 


d^ki d?k2 
(27r)® (27r)3 


5(k{kr)5(k{k2) {0[-kr)0{-k2)) 


+ / n - 3^(0(-fci)0(-fc2)0(-fc3)0(-fc4)) 

+ ^(0(-fci)0(-fc2)Ti,(fci + k2))' {0{-ks)0{-k4.)fkliks + k^))' 
{2^f6^ ( kj^ + k2) }) . 


(5.24) 


From eq.(5.22) and eq.(5.24), we see that the four point function has two types of contri¬ 
butions, 

{5(f)5(f)5(f)54>) = {64>d(j)5(l)6(j))cF + {6cl)5(j)5(j)6(j))ET- (5.25) 


®In [4], the term (TkiOO) in the wave function appeared with a coefficient 1/4 (see eq(2.36) of [4]). But 
in the present work we choose to have a 1/2, which means we need to consistently replace 

{TklOO)there 2{TklOO)here 


while using expressions from [4]. 




















Here, {6(j)5(p64>d(p)cF is the term proportional to (OOOO) 


{6(l){k:i_)S4>{k2)S4>{k3)S4>{k4))cF = 


1 


8 nLi 

and {S(p6(l)S4>6(f>)ET is the term proportional to {OOTij)'{OOT^i)', 

1 


(0(fci)0(fc2)0(fc3)0(fc4)), (5.26) 


{5(t>{kr)5ct^{k2)5(t^{k3)5ct>{k4))ET = 7 


/(fcl,fc2,fc3,fc4) 


+I{k4,k3,k2,k4) +/(fcl,fc4,fc3,fc2) + .^(^3, ^2, ^1, ^ 4 ) 
+ /(fc4,fc2,fc3,fcl) + I{k3,k4,k4,k2) 


(5.27) 


where 7(^1, ^ 2 ; ^ 3 ; ^ 4 ) is given in eq.(E.13) of [4], 

d^ke {fij{k^)fki{kQ)) 


I{k4,k2,k3,k4) = 




(27r)3 (27r)3 /c| A;| (5.28) 

(0(fci)0(fc2)Ti,(fc5))(0(fc3)0(fc4)rfcKfc6)). 


Now, the term {6(f>6(j)64>6(f>)cF is invariant under a special conformal transformation, 
whereas the term {5(p64>5(f)6(j)) et does change. We therefore have 


E2 


a=l 


{S4,{ki)i4,{k2)i<Pik3)S<p{kt)) = 


E4‘. 


a=l 


{6(j){k4)5(p{k2)6(p{k3)64>{k4))ET- 

(5.29) 


As discussed in appendix (E.2) of [4], we have 


E^ 

a=l 


{d4>{ki)6(l){k2)6(l){k3)d(j){k4))ET = 7 ^ 


d^I{k4,k2,k3,k4) 


+(5‘^/(fci, fcs, k2, ^ 4 ) + 5^1{k4,k4, ks, ^ 2 ) + 5^1{k3, ^ 2 , fci, ^ 4 ) 
+ d^I{k4, k2, fcs, fci) + 6^1{k3, k4, fci, ^ 2 ) 


(5.30) 


with (5^J(fci, fc 2 ) ^ 3 ) ^ 4 ) given in eq.(E.23) of [4], 

5^1{k4,k2, ks, k4) = I2h^ I (2^)3 % {0{k4)0{k2)%j{k)) 


(5.31) 


(0(fc3)0(fc4)rfcz(-fc)). 


By using the Ward identity eq.(3.8) of [4] expressed in momentum space, 

k,{f,,{k)o{k{)o{k'^)) = i(A:^,(o(fc; + fc)o(fc')) + fc;,(o(fc'+fc)o(fc;))), ( 5 . 32 ) 

we can calculate the RHS of eq.(5.30). This gives the result eq.(5.21) for the change in the 
four point function, and completes the check. 
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6 Late Time Behaviour of Modes 


In this section, we elaborate on the late time behaviour of modes in the canonical model of 
slow roll inflation and also after including higher derivative terms in the action. 

6.1 The Canonical Model of Slow Roll Inflation 

We have discussed in section 2 that one can use the residual time reparametrization invari¬ 
ance, eq.(2.10), in the gauge eq.(2.2), to set 5(j) = 0, and that the remaining perturbations all 
become time independent in this gauge, at late times. Here we demonstrate this behaviour 
explicitly in the canonical slow roll model of inflation. The behaviour in the presence of 
higher derivative terms is discussed in section 6.2. 

The action for the canonical model of slow roll inflation is given by 




In the canonical model, the Hubble parameter eq.(2.4) is given by 


9 V 
= -. 


( 6 . 1 ) 


( 6 . 2 ) 


The background (/>(t) satisfies the equation of motion 

I + + V'{^) = 0, 

which in the slow roll approximation reduces to 


(6.3) 


(6.4) 


where a ' denotes a derivative with respect to the scalar field. Using eq.(6.2) and eq.(6.4), the 
slow roll parameters ei,(5 and e, defined in eq.(3.4), eq.(3.5) and eq.(3.6), can be expressed 
as 


1 /V 


f\ 2 


u" 


ei = e = - T7 ,<5 = 61 - — 


The slow roll conditions, eq.(3.8), are 


/\ 2 


V 


V" 

« 1 , ^« 1 . 


V 


(6.5) 


( 6 . 6 ) 


For the purpose of convenience in calculations, it is helpful to further decompose the 
metric perturbation jij, eq.(2.5), as follows [41] 


Jij — + didjB + diCj + djCi + Dij], 


(6.7) 


where A, B transform as scalars, Ci transforms like a 3-vector and Dij transforms as a 
rank-2 tensor under spatial rotations. Note that the perturbation Q is divergence-less, and 


26 - 


Dij is transverse and traceless, 


diCi = 0 , diDij = 0 , Du = 0 . ( 6 . 8 ) 

The Einstein equations to linear order in the perturbations about the FRW inflationary 
background are given by (see appendix C)^ 

- '''<« ^ - K^) 

Q = A-a^B -‘iaaB, (6.10) 

— ^S(j) = A, ( 6 . 11 ) 

- (2 <^ 4 - VA^)6(P^ =^A + 3(^^^A + ^V^B+(^^^V^B, (6.12) 

for the scalar perturbations A,B,6(p. The vector perturbations Cj satisfy the equation 

= 0. (6.13) 

For the tensor perturbations Dij we get 

Dij + 3 ~ ~ (6-14) 

The equation of motion for 5(f) is 

6^ + 3(-]6^ + V'U) H = --H^A + V^b) . (6.15) 

\a J 2 V / 

Eq.(6.10) can be used to solve for B in terms of A, 

■ 3 ^Gi(®) + J dt” a{t")A{t", x'^ + G 2 {x), (6.16) 

where Gi,G 2 are arbitrary functions of x. 

The late time behaviour of these equations can be obtained by dropping all spatial 
derivatives of the form V^/a^ in eqs.(6.9), (6.14) and (6.15). In addition, due to the 1/a? 
pre-factor in eq.(6.16), we get that 

B{t,x) ^ G 2 {x) for t —)• 00 , (6.17) 

so that all time derivatives of B vanish at late times. Equations (6.9), (6.12) and (6.14) 
then simplify to 

y'[f/))S(f) = ]^A + 3{^ A, (6.18) 

^For brevity, we present the equations in units with Mpi = = 1. 


B{t, x) = J dt' 


a{t') 
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(6.19) 


- V'{^)dcj)] =-A + 3(-] A, 


Dij + 3 ( — ) Dij — 0, 


and eq.(6.15) becomes 


3 - ; 


5(j) + 3 [ - ] 5(p + V {4>) 64> + -(pA = 0. 


( 6 . 20 ) 


( 6 . 21 ) 


As discussed in appendix C, the late time behaviour for A, 64> is 

A{t, x) = Pi{x) - 2 P 2 ix), 

and 


( 6 . 22 ) 


6(j){t,x) = - 4>{t) P 2 {x), (6.23) 

where Pi,P 2 are time independent functions of x. Also, the perturbations Ci,Dij become 
time independent. 


We can now carry out a time reparametrization 


t^t + P2{x), (6.24) 

along with the accompanying spatial reparametrization, eq.(2.11), which maintains the 
gauge choice eq.(2.2). Note that under the time reparametrization eq.(2.10), and the ac¬ 
companying spatial reparametrization eq.(2.11), the perturbations change as 


M = 2(^) t(x). 

6Ci = 0 , 

5Dij = 0 , 

6{6(j)) = ^e{x). 


(6.25) 

(6.26) 

(6.27) 

(6.28) 
(6.29) 


We see from eq.(6.23) and eq.(6.29) that the change eq.(6.24) sets the late time value of 6(j) 
to vanish. In addition, using eq.(6.25) we see that the value of A is given by 


A^A' = PAx), 


(6.30) 


while C, Dij are unchanged and therefore continue to be time independent. B is changed 
by the the time reparametrization eq.(6.24), see eq.(6.26), but this change vanishes at late 
times, and thus B too continues to be time independent. Thus, we see that in the gauge 
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64> = 0 all the perturbations freeze out at late times 


6.2 Higher Derivative Corrections 

In the discussion above, we have considered the canonical model of slow roll inflation, with 
the action in eq.(6.1). The action for this model involves two-derivative terms. One of 
the main motivations of our work is to be able to use symmetry considerations in more 
complicated situations where explicit computations or models may be unavailable. An 
example is the possibility that the Hubble scale H during inflation is of order the string 
scale Mst, so that higher derivative corrections to eq.(6.1) would be important. Given 
our limited knowledge of string theory in time dependent situations, explicit models or 
calculations for such a scenario are not possible today. But a symmetry based analysis 
should still be possible, as we discuss further in this section. 

The more general situation we have in mind is the one with an effective action having 
higher order terms of the schematic form 


S 


1 


WttG 


j 


R'^ 


(dcj)Y 


R + {d(j)) + + + 


A2 


+ 


(6.31) 


The higher derivative terms are important because H ~ 0(A). In eq.(6.31), A is the 
underlying cutoff scale, which would be of order the string scale Mgt in string theory. The 
term R?/h? schematically denotes four derivative terms, and so on. Also, the coefficients 
of each of the higher order terms can in general be a function of (j). Let us note that in 
the background solution the contribution from terms like will be small, since the 

inflaton will be evolving slowly. However, these terms will be important in determining 
the behaviour of the perturbations, since the perturbations will start out with physical 
wavelengths A <C and then freeze out at a time when A ~ 

Let us note that H ~ 0(A) is consistent with the bounds on the tensor perturbations, 
since A can be much smaller than Mpi, as indeed happens in weakly coupled string theory. 
The condition A <C Mp; also ensures that all quantum loop effects are small, and it is only 
tree level effects involving the higher derivative corrections which are important in the kind 
of scenario we have in mind. 

In fact, considerations of the last few sections can be extended in a straightforward 
way to situations of this type. The crucial point is that even with the higher derivative 
terms present, one can argue that solutions with the same asymptotic behaviour as in the 
two-derivative case continue to exist. The underlying reason for this is that the asymptotic 
behaviour in the two-derivative case follows from gauge invariance. We will discuss this 
in more detail in the next subsection. Given this fact, the arguments leading to the Ward 
identities can be easily seen to apply in cases with higher derivative corrections as well. A 
further change of coordinates allows us to set 5(p to vanish, as discussed in section 2, and 

®The perturbations which appear in section 2 and the discussion thereafter, are given by 

C = id + , and % = -t diCj -b djCi + didjB - i 
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the invariance of the wave function under the residual spatial reparametrizations in the 
synchronous gauge then leads to the Ward identities of interest. 

6.2.1 Freezing of the Perturbations 

We start with a discussion of the spin-2 component Dij, eq.(6.7), which corresponds to grav¬ 
ity waves. In the two-derivative theory it satishes the eq.(6.14). At late times, when k'^ja^ 
becomes sufficiently small, this becomes eq.(6.20), which has the general solution eq.(C.18). 
In particular, Dij becomes time independent, satisfying eq.(C.19), since the additional so¬ 
lution proportional to Kij in eq.(C.18) dies out as t —)■ oo. Higher derivative terms would 
result in contributions to the equation of motion with either additional spatial derivatives, 
and/or additional time derivatives. All terms with spatial derivatives will become small, 
since the physical wavelength A for fixed k becomes large at late times. Thus the only 
terms which survive will have additional time derivatives. It is then clear that the solution 
eq.(C.19) will continue to hold even when higher derivative corrections are included. 

However, in the presence of higher derivative terms there could be additional solutions 
which do not die out at large t. We will assume that the correct boundary conditions in 
the far past are such that any such solution is not “turned on” in the far future, leading to 
eq.(C.19). 

Exponentially growing solutions would signify an instability. Our assumption that 
they are absent is consistent with the background inflationary solution being stable. There 
could be additional oscillatory solutions though, which are non-decaying. We cannot rule 
these out except by appealing to the initial conditions. However, the following possibility 
is worth mentioning in this context. The additional oscillatory solutions might be present 
if the higher derivative corrections in eq.(6.31) arise in the first place by integrating out 
massive particles with a mass ~ 0(A). This could happen in an underlying theory where 
all particles, the massive ones and the graviton, satisfy second order equations of motion, 
leading to a well posed initial value problem. In this case, the graviton will indeed have 
the solution discussed above, eq.(C.19), with a second solution which decays, eq.(C.18). If 
A ~ 77, these additional particles would also be produced during inflation, with a suitable 
Boltzmann suppression. However, the Ward identities we derive in section 3 will continue 
to hold in this case as well. The wave function in the presence of these fields will continue to 
be invariant under spatial reparametrizations, and thus after integrating these heavy fields 
out, eq.(3.9), the same Ward identities will follow for ^ and jij. 

The discussion for spin-1 is even more straightforward. The solution found in the two- 
derivative case is pure gauge, since there are no physical degrees of freedom with spin 1. 
Starting from the unperturbed solution of the form eq.(2.3), and carrying out a spatial 
reparametrization 

X* —)• x*-|-e*(®), (6.32) 

one gets 

Ci = Ei - did~‘^{d ■ e), (6.33) 

so that the most general time independent Ci can be turned on with a suitable choice of 
e*(a;). This makes it clear that a solution of the form eq.(C.17) must continue to exist in 
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the presence of higher derivative terms too. 

Finally we come to the scalar perturbations. In the two-derivative theory, the late time 
behaviour for solutions was found to be eq.(6.22), eq.(6.17) and eq.(6.23), for A,B and 6(j) 
respectively. We now argue, in analogy with the case of spin-1 above, that the existence of 
solutions exhibiting this behaviour follows from spatial and time reparametrizations which 
preserve the synchronous gauge eq.(2.2). Starting from eq.(2.3) and doing the transforma¬ 
tion eq.(6.32) gives 

B = 2d-^{d-e), (6.34) 

so that the most general time independent B can be turned on. Also, starting from eq.(2.3) 
and carrying out a transformation 


x^-^x^{l + e), (6.35) 

where e is a constant, gives 

A = 2€. (6.36) 

The late time behaviour of A with higher derivative terms will still be determined by an 
equation where all spatial derivatives can be dropped. The solution eq.(6.36) then means 
that actually 

A^Pi{x) (6.37) 

will be a solution to the small perturbation equations. 

Finally, doing the time reparametrization eq.(2.10) gives rise to the solution 

A = 2He{x), 64> = ^€{x). (6.38) 

Putting all these solutions together, we get the general late time behaviour seen in eq.(6.22), 
eq.(6.17) and eq.(6.23). 

Since these solutions in the spin-0 case arise just from gauge invariance, they will 
continue to be true even in the presence of the higher derivative terms. As in the case of 
the spin-2 mode, there could as well be additional solutions which do not decay, but we will 
assume that they are either not turned on due to the initial conditions, or are of oscillatory 
type arising due to additional massive particles, which do not invalidate the arguments for 
the Ward identities. 

7 Conclusions 

In this paper, we have derived the Ward identities that arise from scale and special conformal 
transformations for single held inhation. Our results are given in eq.(3.25) and eq.(3.38) for 
the scalar perturbations, and eq.(3.27) and eq.(3.41) for the tensor perturbations. Similar 
results for mixed correlators can also be easily obtained, see eq.(3.28) and eq.(3.44). 

The Ward identities for the special conformal transformations also involve a contri¬ 
bution due to a compensating spatial reparametrization, as explained in section 3.2. The 
underlying reason for this is that we are working with local correlators in a quantum theory 
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of gravity. Such correlators can be defined in perturbation theory after suitable gauge fix¬ 
ing, but a compensating spatial reparametrization must then be carried out to preserve the 
gauge, for deriving the Ward identities of special conformal transformations [4]. The Ward 
identities for scale invariance do not require such a compensating transformation and are 
well known in the literature already, [1], and called the Maldacena consistency conditions. 

The Ward identities we obtain also incorporate the breaking of the symmetry. 

In fact, this breaking is incorporated to all orders in the slow roll parameters. The resulting 
relations can be thought of as being the analogues of the Callan-Symanzik equation, but 
now for both scale and special conformal transformations. 

When the slow roll conditions are approximately valid, the Ward identities impose use¬ 
ful constraints on the correlation functions. The coefficient functions which appear in the 
wave function, and which transform in a manner analogous to correlation functions in a 
conformal field theory, can then be constrained order by order in the slow approximation, 
and the resulting constraints on the expectation values, in agreement with the Ward iden¬ 
tities, can then be obtained. For the scalar three point function, which is observationally 
the most important one for non-Gaussianity, this was discussed in [5]. 

We work in a theory where the degrees of freedom are the metric and a scalar field. 
However, it is worth emphasizing that our results are also valid in situations where there 
are extra massive fields present during inflation, with masses of order the Hubble scale, or 
even higher. The wave function in the presence of such fields must still meet the equations 
of motion imposed by varying the shift and lapse functions, and thus must be invariant 
under the spatial reparametrizations discussed in section 4.1, in the gauge eq.(2.17) at late 
times. As a result, after the heavy fields are integrated out in deriving the expectation 
values for ^ and %j, in the step analogous to eq.(3.9), the same Ward identities as before, 
eq.(3.25), eq.(3.27), eq.(3.38) and eq.(3.41) are obtained. 

In contrast, our results are not valid when there are additional scalar fields which are 
much lighter than the Hubble scale, as in multi-held models of inflation. In this case, it is 
well known that the results are model dependent, and do not follow just from the underlying 
symmetries. 

It will be worth examining the Ward identities of scale and special conformal transfor¬ 
mations also in situations where the conformal symmetries are badly broken. This happens, 
for instance, in DBI inhation [42], [43]. It can also happen if the initial state is not the 
Bunch-Davies vacuum, and breaks the conformal symmetries to a signihcant extent. 

In the context of AdS physics, space-times where conformal symmetry is broken are 
important in the study of condensed matter physics and QCD. Examples include Lifshitz 
and hyperscaling violating geometries. The Ward identities for the stress tensor etc. can 
be obtained in such situations in a way completely analogous to what we have used above. 
Some discussion of the identities in such situations can be found in [44]. 

The Scalar Three-Point Function '. Since the scalar three point function is of the great¬ 
est interest as a test of non-Gaussianity, let us end by commenting on it in some more detail. 
The Ward identities of interest here are eq.(3.25) and eq.(3.38) for n = 3, and relate the 
scalar 3 and 4 point correlators. These Ward identities were studied to the leading order in 
the slow roll expansion in [5]. The resulting relations in terms of coefficient functions are 
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given in eq.(3.24) and eq.(3.25) of [5] 


a=l 


d 

dka 


(0(fcl)0(fc2)0(fc3)) = ^ {0{ki)0{k2)0{ks)0{k4)) 


(7.1) 




and 


cU0{kMk2)0{ks)y + cUO{k^)0{k2)Oiks)y + £L(0(fci)0(fc2)0(fc3))' 


= 2 - 


H 


b- 


d 

dka 


{0{k^)0{k2)0{ks)0{k^)y 


fc4—>-0 


(7.2) 


with defined in eq.(4.22). Note that in the leading slow roll approximation, the four 
point coefficient function (OOOO) can be calculated in the conformally invariant limit [4]. 
As a result of the factor of cj) on the RHS of eq.(7.1) and eq.(7.2), the three point function 
{OOO) will be suppressed. Converting to expectation values, one gets that 


(CCO 

(CC)2 


(7.3) 


where the slow roll parameter e is given in eq.(3.6). Although the functional form one will 
get in general is different, this roughly corresponds to 


Inl ~ e- 


(7.4) 


ft is well known that the parameter r which measures the ratio of the power in the 
tensor to scalar perturbations is given by ® 


Pcik) 


16e. 


(7.5) 


Note that in theories which are not of the type described by a canonical model of 
inflation, eq.(6.1), e.g., those involving higher derivative corrections, ei and e as defined 
in eq.(3.4) and eq.(3.6) need not be the same, fn these theories also, eq.(7.3) with the 
definition of e given in eq.(3.6) is still valid to leading order in the slow roll parameters. We 
see from eq.(7.4) and eq.(7.5) that there is therefore an interesting tie-in between the ratio 
of power in the scalar and tensor perturbations, and the non-Gaussianity. This connection 
is well known in the canonical slow roll models, but we see here that it is more general, 
since eq.(7.4) follows from symmetry considerations alone. 

The estimate in eq.(7.3) should actually be thought of as a lower bound. A contribution 
due to an intermediate graviton (or the stress energy tensor running as an intermediate 
in the {OOOO) correlator) will give a contribution of this order to the non-Gaussianity. 


®The tensor and scalar power spectra Pt(fc),Pc(fc) are 


Pt{k) 


4 


and Pcik) 


1 1 
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However, as has been emphasized in [45], if there are additional particles of mass of order 
the Hubble scale which couple more strongly than the graviton, the contribution can be 
even bigger 

Keeping the above considerations in mind we can phrase this tie-in between the two 
scales as follows. If tensor perturbations are observed in the future, so that e is known, 
we would have a firm prediction on a lower bound on non-Gaussianity that follows only 
from conformal invariance. On the other hand, if the non-Gaussianity is observed and 
found to be of a bigger magnitude than the bound on e that arises from constraints on the 
tensor perturbations, eq.(7.5), then it would rule out the scenario of approximate conformal 
invariance. More correctly, it will rule out this scenario together with the assumption that 
particles which appear as intermediate states, and contribute to the non-Gaussianity, couple 
to the inflaton only with gravitational strength. 
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A Transformation of Perturbations under Spatial Reparametrizations 

In this appendix, we would like to give some details about the transformation properties of 
the perturbations under spatial reparametrizations. We consider the perturbed line element 
in the gauge eq.(2.2), 

ds'^ = — dt^ + hij{t, x) dx^dx^, (A.l) 

with 

hij = a^{t) Qij = a^{t) [6ij + %•], (A.2) 

where 

% = 0. (A.3) 

Consider now a spatial reparametrization of the form eq.(2.7). The change in hij under this 
transformation is 

5hij = ViVj + VjVi. (A.4) 

^^Similarly, in theories where conformal invariance is violated to a significant extent, the non-Gaussianity 
can be bigger, e.g. in DBI inflation. 
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Eq.(A.4) implies 


^9ij 


1 

1 


[diVj + djVi - 2 r“- 

[hjkv^^ + {hikv^^ 


gjk div^ + Qik djv^ + v^dkQij , 


-Va h‘^^{dihji, + djhii, 


dbhij) 


where indices will now be raised and lowered by 5ij. Eq.(A.5) gives us 

^9ii — ‘^9ikdiV + V dk9ii- 


(A.5) 


(A.6) 


Putting Qij from eq.(A.2) in eq.(A.6) gives the change in under spatial reparametrizations, 
eq.(2.7), as 

^ dm + v^dkC + ^ diVj , (A.7) 

which is the result quoted in eq.(4.1). Once we have calculated we can insert the full 
Qij in eq.(A.5) to get the change in qjj as 


= { diVj + djVi - - daVa Sij ) + { 7ifc djV + ^jk diV + 


2 2 \ 

+ V dk'Jij ~ daVa ^ij ~ daVf) ^ah i^ij 4“ 'Jij') j ■ 


(A.8) 


For simplicity, we call the terms in eqs.(A.7) and (A.8) which are proportional to the 
perturbations as the homogeneous pieces of the transformation, and the parts independent 
of the perturbations as the inhomogeneous pieces of the transformation. 

Having obtained eq.(A.7) and eq.(A.8), we can calculate the changes SC and Sjij for 
the specific cases of scale transformations, eq.(3.12), special conformal transformations, 
eq.(3.29), and the compensating spatial reparametrization, eq.(3.33). For scale transforma¬ 
tions, the change in C and is given by eq.(3.13) and eq.(3.14) respectively. Similarly, for 
the special conformal transformations, the changes are given by eq.(3.30) and eq.(3.31), and 
for the compensating spatial reparametrization, the changes are eq.(3.34) and eq.(3.35). 


B The Scalar and Tensor Spectral Tilts 


As a simple check on the Ward identities, we consider here the 2-point correlators. For 
scalar perturbations, the scaling Ward identity in eq.(3.25) relates the 2-point expectation 
value to the 3-point expectation value. 


d 


3+E j = - (c(mcFw 


(B.l) 


fes—>^0 
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The expression for (C(^i)C(^ 2 )) is 


1 


(C(fci)C(fc2)) = {27:f5\kr + ^ ^ 


(B.2) 


where ns is the scalar tilt. Thus, we get 
2 


.a=l 


d 

dka 


(C(fci)C(fc 2 ))' = (-3 + ns)(C(fci)C(-fci))', 


(B.3) 


which on substituting back into the eq.(B.l) gives the well known Maldacena consistency 
condition 


lim (C(fci)C(fc2)C(fc3))' = -ns (C(fci)C(-fci))' iCiksKi-ks))'. 


(B.4) 


By using the expression for (('(fci)^(fc2)C(^3)) from [1], we get 
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lim {C{ki)C{k2)C{k3)y = (6e - 2ry) (C(A;i)C(-fci))' (Cik^Xi-ks))'. (B.5) 

fcs-^O 

Putting the limit from eq.(B.5) into eq.(B.4) gives the expression for the scalar tilt as 

ns = 27] — 6e , (B.6) 

which is indeed the correct expression, [1]. 

Similarly, consider the tensor Ward identity in eq.(3.27), with n = 2. This has the form 




. (B.7) 


fcs—>-0 


In writing eq.(B.7), we have introduced the two polarization tensors for the graviton, ef^, 
through the relation 


%j{k) = '^elj{k)^s{k). 


(B.8) 


S = 1 


Now, the expression for ( 7 s(fci) 7 s'(fc 2 )) has the form 




{%iki)7s'ik2)) = {27rf 6^{k^ + fca) 5.,.' ^ 
where nx is the tensor tilt. By using the expression eq.(B.9) in eq.(B.7), we get 


(B.9) 


lim (7s(fci)76'(^2)C(fc3))' = - nr 5s,s' (7s(^i)7s(-fci))'(C(fc3)C(-fc3))'- (B.IO) 
We can calculate the limit on the left side of eq.(B.lO) by using the expression for the 


^The slow roll parameter 77 is given by 77 = e — 5 = 
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correlator (7s(fci)7s'(fc2)C(^3)) from [1], This gives 


lim {%iki)%>{k2)Cik3)y = ‘2e6s,s'{%iki)%i-k-i_)y{C{k3)({-k3)y. (B.ll) 

Then by comparing eq.(B.lO) and eq.(B.ll), we get 

nT = -2e, (B.12) 

which is the correct expression for the tensor tilt. 


C The Behaviour of Perturbations in Canonical Slow Roll 

In this appendix, we provide some details of the analysis given in section 6.1. We follow [41] 
for our calculations. Our gauge choice, eq.(2.2), is same as the synchronous gauge of [41] (see 
section 5.3 (B)). The relevant equations are eqs.(5.3.28)-(5.3.33) for scalar perturbations, 
eq.(5.1.51) for vector perturbations, and eq.(5.1.53) for tensor perturbations. 

The energy-momentum tensor for the inflaton can be calculated by varying the matter 
part of the action (6.1) with respect to the metric. It is given by 

^ (C.l) 

This has the form of the energy-momentum tensor for a perfect fluid, 

= {p + P) u^u'^ + Pg^'' , (C.2) 

with the energy density p, pressure P, and the four-velocity given by 


p = -\iv<py + vicp), 

(C.3) 

P=-^(V0 )2-P(^), 

(C.4) 


(C.5) 


For our purpose, we specialize to the case of single field slow roll inflation. We then 

have 


p= + Vi(p), 

(C.6) 

P = ^$^-V{^), 

(C.7) 

yP = !,«* = 0, 

(C.8) 


for the homogeneous background (p{t). By expanding eqs.(C.3)-(C.5) to linear order in the 
perturbation 5(j), we get 

5p = ^>6(1) + V'{4>) '5'^) 


(C.9) 


SP = 


du = - 


(C.IO) 

(C.ll) 


where 6u is defined through 6ui = didu + SuY, and 6uY = 0 for single held inhation. Also, 
for single held inhation, the anisotropic stresses in the perturbed energy-momentum tensor 
vanish, 

(C.12) 


= 0, ttY = 0, vrh = 0. 

By using the eqs.(C.6)-(C.12) above in the eqs.(5.3.28)-(5.3.33), eq.(5.1.51) and eq.(5.1.53) 
of [41], we obtain the perturbed Einstein equations eqs.(6.9)-(6.15) given in section 6.1 for 
the scalar, vector and tensor perturbations, along with the equation of motion for the back¬ 
ground eq.(6.3). Note that the perturbations Gj in eq.(5.1.51) of [41] vanish due to 
our gauge choice eq.(2.2). 

We now provide some details for calculating the late time behaviour of the perturba¬ 
tions. To solve for the perturbation A, we consider eq.(6.18). Inserting 6(j) from eq.(6.11) 
into eq.(6.18), we get an equation purely for the perturbation A, 


A 

2 ^ 


31 “') + '"'® 

a 


A = 0. 


(C.13) 


By using the background eq.(6.3) in eq.(C.13), we get 


A-2 


A = 0. 


(C.14) 


The general solution to eq.(C.14) is 

A{t,x) = Pi{x) + P 2 {x) j 


(C.15) 


where Pi{x), P 2 {x) are two arbitrary functions of x. Eq.(C.15) on using the background 
equation 


H = 


d f d 


dt \ a 


1 


(C.16) 


becomes 


A{t, x) = Pi{x) - 2{- ] P 2 {x), 


which is the solution quoted in eq.(6.22). 

Once we have obtained the solution for A, it is straight forward to obtain the solution 
for the perturbation 64>. From eq.(6.11) and eq.(6.22), it follows that 

5(t){t,x) = -${t)P2{x), 

as given in eq.(6.23). One can check explicitly that the solutions eq.(6.22), eq.(6.23) satisfy 
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the other equations, namely eq.(6.19) and eq.(6.21). 

The equation for the perturbation Cj, eq.(6.13), has the general solution 


Ci{t,x) = d ^Qi{x) 


(C.17) 


which shows that the perturbation Ci is frozen for non-zero momentum modes, which are 
the ones of interest to us. 

Finally, we consider eq.(6.20) for the tensor perturbations. The general solution is 





which in the late time limit also gets frozen. 


Dij{t, x) K, Dijix) for t —>■ oo. 


(C.18) 


(C.19) 
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